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We report Keck High Resolution Echelle Spectrometer data and model atmosphere analysis of two 
helium-dominated white dwarfs, PG1225-079 and HS2253+8023, whose heavy pollutions most likely 
derive from the accretion of terrestrial-type planet(esimal)s. For each system, the minimum accreted 
mass is ~10 22 g, that of a large asteroid. In PG1225-079, Mg, Cr, Mn, Fe and Ni have abundance 
ratios similar to bulk Earth values, while we measure four refractory elements, Ca, Sc, Ti and V, all at 
a factor of ~2-3 higher abundance than in the bulk Earth. For HS2253+8023 the swallowed material 
was compositionally similar to bulk Earth in being more than 85% by mass in the major element species, 
O, Mg, Si, and Fe, and with abundances in the distinctive proportions of mineral oxides - compelling 
evidence for an origin in a rocky parent body. Including previous studies we now know of four heavily 
polluted white dwarfs where the measured oxygen and hydrogen are consistent with the view that the 
parents bodies formed with little ice, interior to any snow-line in their nebular environments. The growing 
handful of polluted white dwarf systems with comprehensive abundance measurements form a baseline 
for characterizing rocky exoplanet compositions that can be compared with bulk Earth. 
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G 241-6) - Stars: abundances - white dwarfs 



Stars: individual: (HS2253+8023, PG1225-079, 



INTRODUCTION 



Accretion by white dwarfs of tidally disrupted bodies from their planetary systems ([Jural 120031 ) has become 



the most viable explanation for the presence of high-Z (atomic number Z>2) elements i n the vast majority of 



olluted PAZ, DB Z, and DZ white dwarf atmospheres with Teff < 20,000 K ([Jura 



2008 



Zuckerman et al 



2010; 



arihi et al.ll2010al . and references therein). The basic idea is that smaller planetary bodies such as asteroids or dwarf 



^planets, can acquire highly eccentric orbits by perturbations from larger planetary bodie s, especially considering th e 
dynamical re a rrang ement of the planetary system during stellar evolution and mass loss (jDebes fe Sigurdssonl l2002). 



Bonsor et al.l (|201ll ) demonstrate that the post-stellar mass loss evolution of a planetesimal belt can supply the 
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required mass influx to explain observed polluted white dwarf (WD) accretion rates. Objec ts that .jou r ney w ithin 
the Roche radius of the WD will experience tidal shredding, forming a disk of dust and/or gas (j Jurall2003 . 2008 ). and 
subsequently accrete into the star's, otherwise pure, hydrogen and/or helium atmosphere. This natural mechanism 
offers the potential to measure the distilled elemental constituents of a planetary parent body or bodies - a unique 
and powerful tool, indeed. 



High-resolution and high-sens i tivity spectroscopy, along with state-of- t he-art model atmospheres (|Koesterll2009i . 



20101 ; iTremblav Bergeron 



2009 



Dufour et al 



2010 



Vennes et al. 



2010, 



2011 



), has opened the door to detailed 
studies of WD atmospheres and the high-Z material that pollutes them. To date, a handful of WDs with very 
rich spectra have been studied at high-resolution and high-sensitivity, beginning with the Keck/ HIRES follow-up 



Zuckerman 



et al.l 120071) of GD 362 , whose extreme high-Z-rich atmosphere was pointed out by iGianninas et al 



mm and iKawka Yenned (120061 ). The d etection of oxyge n (O) along with the other major elements (Mg, Si, 
and Fe) in GD40's polluting parent body (jKlein et al.l |2010| ) established chemical evidence - the stoichiometric 



balance of mineral oxides 



(Vennes et al. 


2010: 


Dufour et al. 


2010; 



that externally polluted white dwa rfs have ent e red a stage of being practical measuring devices of the elemental 



compositions of rocky exoplanets. iBond et al.l (|2010l ) predict a wide diversity of terrestrial exoplanets including 



Earth-like, refractory-rich, and carbon-rich obj ects. Polluted white d warfs provide real data at the level of detail 



needed to make discerning measurements (e.g.. IZuckerman et al.ll2011 



In this paper we report the detailed analysis of HIRES spectra of HS2253+8023 (V mas =16.1, DBAZ, T e ff ~ 
14,400 K; hereafter HS2253) and PG1225-079 (=K 789-37, V mas =14.8, DZAB, T e// ~ 10,800 K; hereafter PG1225). 
Both stars were previo usly known to displ a y features of H , Mg, Si, Ca and Fe ma inly from low (~ 6A) resolution 
IUE and optical data (jKoester et al.lll99ol ; iFriedrich et aDll999l ; IWolff et al.l hooj). Ca lcium in PG1225 was also 
(re)observed and analyzed from the ESO SN la progenitor (SPY) survey (jKoester et al.ll2005bl ). The HIRES obser- 
vations provide a new and important detection of oxygen in HS2253, completing the set of major terrestrial elements 
- O, Mg, Si, and Fe - as well as H, Ca, Ti, Cr, Mn and possibly, Ni. For PG1225 the HIRES data display six newly 
identified elements: Sc, Ti, V, Cr, Mn, and Ni. By analyzing all the detected elements together in a given model, 
with multiple lines for each element, and accounting for abundance sensitivity to errors of the model parameters, 
we calculate pollutant-to-pollutant abundance ratios with an estimated uncertainty of ~ 40% or less, a significant 
improvement over previous analyses in which derived abundances relative t o helium are typica lly uncertain by 0.2-0.7 



dex. Infrared studies show that HS2253 does not have an infrared excess (Farihi et al 



a subtle infrared excess appearing only at 7.9/im (jFarihi et al.ll2010bl : iKilic et al.ll2008al ). Notwithstanding the lack 



20091 ). and PG1225 has just 



of detected circumstellar dust, we demonstrate that these systems have accreted an appreciable amount of material 
with terrestrial planet (esimal) composition. 

We also report the detection of Ca II HfeK emission cores in PG1225, the first time seen in a white dwarf. 
Other peculiar line profiles observed with HIRES are narrow absorption cores in helium lines of PG1225, GD362 
and GD16, which are displayed in this paper. 
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There are a number of benefits of working with helium-dominated WDs. One is that, compared to hydrogen- 
dominated WD atmospheres, the lower opacity helium environments more readily display tiny concentrations of 
numerous contaminant elements, e.g. the trace elements vanadium and scandium are detected in PG1225 at 10 - 
11 orders of magnitude lower abundance than the helium background they are diluted in. Secondly, det ailed model 



atmo s phere analysis of he lium-atmosphere WDs indicates the hydrogen in these stars must be accreted ()Voss et al 



20071 : iDufour et al.l 120071 ); thus when helium is the background constituent, one obtains a measure ment of the 



polluting hydrogen abu ndance, which is important for understanding the potential accretion of water ice ([Jura et al 



2009 ; I Jura &: Xull2010i ). A simplifying aspect of the helium-atmosphere WDs is that in our T e jj range of interest (< 
20,000 K), all have outer convection zones where the elements are homogeneously mixed. Finally, another benefit 
of working with helium-rich WDs is that with relatively large convective zones (relatively long settling times), the 
pollutants are not disa ppearing so fast from our view. Thus, a higher fraction of polluted WDs can b e found (e.g 



Zuckerman et al 



(2010) 



and 



2010), a nd me aningful lower limits on accreted masses can be obtained. Recently, IDufour et al. 



Koester et al.l (120111 ) have established that absolute minimum accreted masses for a number of polluted 



WDs are >10 23 g, close to that of Ceres, the largest asteroid in our own Solar system. 

Gravitational settling times are always less than the age of the WD, but i n helium-atmo sphere stars of ^eff < 
20,000K, the settling (diffusion) time scales for most elements are 10 3 -10 years ( Koesterll2009l ) . Thus a measurement 
of atmospheric abundances in a helium-dominated WD is a snapshot of one s tage of the evolving interplay of accretion 
and diffusion in that system. As was found for GD40 ([Klein et al.l l2010i ). the patterns of pollutant abundances 
themselves, can narrow the range of possibilities for the accretion-diffusion state of a WD, and we utilize some such 
constraints in the present work. 

The article is laid out as follows: observations, data reduction, and measurements of the spectra are described in 
Section [2 In Section [3] we discuss effective stellar temperature, surface gravity, and the synthetic model atmospheres 
from which we derive the polluting element abundances presented in Section SJ We describe the effects of accretion 
and diffusion in Section [5j incorporate that into our results - including results of previous studies - in Section [61 and 
give our conclusions in Section [71 



OBSERVATIONS & MEASUREMENTS 



The data were acquired with the High Resolution Echelle Spectrometer (HIRES. IVogt et al.lll994l ) on the Keck 
I telescope of Mauna Kea Observatory. Table Q] is the observation log. Typically, a l'/148 slit width was used, 
providing a spectral resolution of R = A/AA ~ 40,000. To compensate for inferior seeing on 2008 Nov 15, a l'/722 
slit was used, resulting in a resolution of R ~ 28, 000 for those spectra. Before combining exposures of the same star 
with differing setups, the higher resolution spectra were degraded by Gaussian convolution to approximate the lower 
resolution spectra. 



Data reduction and absorption line measurements were carried out using PyRAF in a similar manner as described 
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for GD40 in lKlein et al.l (|2010l ). For the present work we took some extra steps, since in addition to a fully continuum- 
normalized version of each spectrum, we also wanted a version that preserved broad-line (> 1-2 A) features. For 
this purpose we used observations of calibration white dwarfs (BD+28 4211, G191-B2B, or Feige 34) in the flux- 
calibration p rocedure to align t he echelle orders (however, absolute flux calibration is not done). As is known for 
HIRES data (ISuzuki et al.ll2003l ). echelle artifacts irregularly appear as bends or waves in the continuum. We visually 
inspected and compared regions of similar wavelength coverage, and selectively combined only those portions of the 
spectra that appeared not to suffer significant distortions. Mostly this amounted to omitting some echelle-order 
edges. Figures Q] and [2] show the high-quality output, and success with broad lines, from these procedures. Residual 
artifacts are low-level and fairly broad, which for example can be seen as waves in the continuum in the region 
A4550-4680A of Figure [2 

For HS2253 the continuum signal-to-noise (SNR) of the data ranges from 25-40 for wavelengths in 3400-8000A, 
and is 10-15 at wavelengths shorter and longer than that. The spectrum of PG1225 is very high quality with the 
SNR of the unsmoothed continuum ranging from 80-110 for wavelengths in 3200-6000A and 40-70 at wavelengths 
shorter and longer than that. 

The line identifica tions and equivale nt widths are listed in Table [21 and were carried out in a similar manner as 
described for GD40 in lKlein et al.l (|2010i ). For PG1225 the many broad lines and high number density of lines was 
not amenable to our usual method of measurement error estimation. Instead, based on our experience with HIRES 
data and checks over multiple exposures, for PG1225 we set the EW measurement errors to 20% for most lines, and 
30% or more for a few lines where the continuum or interval surrounding the line is especially uncertain. We have 
also included here the details of the line measurements of G241-6 (SDSS g = 15.5 mag, DBZ, TW f=1 5,300 K, log 
5=8.0), whose high-Z rich nature and pollutant abundances were reported by IZuckerman et al.l (|2010l ). Especially 
in the case of the very rich spectrum of PG1225, there are a number of blended lines, which are indicated with a 
'b' superscript in Table [2j In a few regions, such as the Mg I A3835 triplet, instead of taking the EW from fitting 
a Lorentzian or Voigt profile, we report the EW from a numerical sum of the flux under the continuum for a given 
interval surrounding the features. Owing to the excellent SNR of these spectra, the measured line strengths span a 
remarkably large range - a factor of 3000 - from 0.011 A (V II), to ~1 A (Ti II) to greater than 30 A (Ca II H&K 
blend). All observed absorption lines, and those not observed, are consistent with expectations of line strengths for 
our derived abundances and SNRs of the spectra. Displays of various portions of the HIRES spectra are given in 
Figures [U-[I3 

For each transition, the shift of line center from the laboratory wavelength is measured, and converted to a 
heliocentric velocity. We use the average line velocity for a systemic value, and Doppler shift the models to the 
measured frame of the WD, which are plotted in red on the figures of spectra. The net systemic heliocentric 
velocities, including the stars' gravitational redshifts, are 4 ± 4 km s" 1 for HS2253 and 49 ± 3 km s" 1 for PG1225. 
From the histogram of line velocity measurements shown in Figure HU for PG1225 we found a range of outliers (filled 
bars in the figure) that represent a detectable shift of the set of Fe I lines (see also Figure fl~2"j) . These are not of 
an interstellar origin, since as indicated in Table [2j most transitions arise from excited atomic levels. An analogous 
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shift for Si II lines has been examined and modeled with Stark shifts bv lVennes et al.l (|201ll ). thus it is conceivable 
that a similar process is responsible for the effect on Fe I in PG1225; further investigation remains for future work. 
HS2253's velocity measurements do not indicate any distinctive sub-populations. 

With such high-sensitivity observations of PG1225 a new feature has turned up (literally) that we have not 
seen before in WDs - emission cores in the Ca II H&K lines as shown in Figure [15j The emission features are 
seen in all 3 separate exposures covering the Ca H&K region of PG1225, with no noticeable variations between 
exposures. Emission cores have been seen in Ha of hydrogen-dominated white dwarfs hotter than 25,000K, and are 
understood to arise from non-local-thermodynamic-equilibrium (NLTE) effects. The one HIRES exposure covering 
Ha in PG1225, taken about 2 weeks later, does not display clear Ha core emission with a continuum SNR~50 near 
6650A. Still, NLTE effects could conceivably cause emission in the strong Ca H&K lines since the cores are formed 
high up in the atmosphere, while the relatively weaker Ha line of PG1225 is likely formed at a deeper location. 
Alternatively, Ca H&K emission cores, as well as Ha core emission, are well known signatures of chromospheric 
activity i n main sequence star s; further discussion of chromospheric activity in a WD is beyond the scope of this 
work (see iMusielak et al.ll2005l ). 



One more spectroscopic peculiarity we point out is the presence of a narrow absorption core component at the 
line center of He I A5876 in PG1225, as well as in GD362 and GD16 (= HS 0146+1847) shown in Figure QB These 
three WDs have helium-dominated at mospheres, with effe c tive temperature s in the range 10,50 K - 11,500 K; they 
all have d etected circumstellar dust (Beckli n et alJ 120051 : iKilic et al.1 120051 : iFarihi et al.l 120091 ) . heavy atmospheric 
pollution (jKoester et al.ll2005al : IZuckerman et al.l 120071 ). and possess relatively large amounts of hydrogen in their 
convective envelopes. Two possibly analogous, but weakly detected, features in PG1225 are narrow absorption lines 
instead of broad ones centered on He I A3889 and He I A 4472, with the latter sh own in Figure [HJ The He I A5876 core 



featur e was first discovered in HIRES data of GD362 bv lZuckerman et al.1 (120071) and noted ag ain bv lTremblav et al 



(|20ld ): we now see it in HIRES spectra of GD16 and PG1225 as well. IZuckerman et al.l (|2007l ) suggested that NLTE 
effects, similar to those seen in the cores of Ha lines in hydrogen-dominated white dwarfs, could be present. The 
exact origin of the He core feature is not known at this time and we do not discuss it further in this paper. 



MODEL ATMOSPHERES 



The model input physics is described by iKoesterl ([20101 ) with all absorption line data taken from the Vienna 
Atomic Line Database (VALD). Key input parameters are the effective temperature, T e jj, surface gravity, g, and the 
atmospheric abundances of elements relative to the dominant constituent, in the current study, helium. In the range 
of temperatures considered here, helium is mostly neutral, and the high-Z elements and hydrogen are significant 
contributors of electrons in the atmosphere. Thus, we calculated and analyzed our models with hydrogen and high-Z 
elements in a self-consistent way, and included Earth-like abundances for undetected but expected elements (such as 
O in PG1225). More than 2000 continuous absorption coefficients were incorporated to account for the blanketing 
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effect. 



In tuning the model inputs to reproduce an observed WD spectrum, it is desirable to achieve consistency for 
Teff, log g and abundance estimates that derive from 1) photometric fitting of the spectral energy distribution 
(SED), preferably using UV to optical and infrared wavelengths, 2) fitting the hydrogen and helium lines in profile 
and/or strength, and 3) demanding a balance in the abundances that come from different ionization states of the 
same element. In principle, the model would meet the requirements of 1) through 3) while also achieving a reasonable 
fit of the broadening and symmetry of heavy element lines - in practice, it can be challenging to reconcile all these 
points, but we find reasonable consistency for the WDs studied here. Precise parameters are not required for our 
abundance analys is, since we are in terested in the pollutant-to-pollutant ratios, which are relatively insensitive to the 
particular model (jKlein et al.ll20ld . also Section [4] of this paper). 



In the analysis of helium dominated atmospheres with effective temperature below ~16,000 K, we must deal with 
the difficulty that the parameters, T e jy and gravity, <?, are effectively coupled, such that T e jj/log g fits to the helium 
lines are not unique. This has led many authors (present authors included) to set log g to the average value of 8.0, 
so that one may proce ed with a fit to the temperature. In the current analysis we took a different approach, similar 



to lDufour et al.1 (120101 ) . who found that a larger than average gravity was required to fit simultaneously the features 
from two ionization states, Mg I and Mg II, in the extremely polluted helium-rich WD, SDSS J073842. 56+183509. 6 
(hereafter, SDSS0738). For each star in the present study we allow for a departure from the average gravity value, 
and instead work with a set of 3 T e jj/log g solutions for T e ff and log g which do a reasonably good job of satisfying 
the spectral energy distribution (SED), He I line strengths, and the ionization balance between neutral and singly 
ionized species of Fe, Mg and Ca. The model parameters used are listed in Table [3] and span a range of ±300 or 
400 K for effective temperature and ±0.3 dex for lo g q. An ad d itional constraint on the gravity could come from 



a sufficiently precise parallax measurement - e.g. as iKilic et al.l (|2008bl ) carried out for GD 362 - and help narrow 
the range of possible log g and T e ff values. To our knowledge, PG1225-079 and HS2253±8023 do not have reported 
parallaxes at this time 3 . 



3.1. HS 2253+8023 T e// /log g 



The SEDs of our current models for HS2253 in the range T e// =14000 - 14 800 K, scaled to 2MA S S pho tometry, 
are in agreement with the IUE data and GALEX photometry (Figure [TTJ) . As iLimoges fe Bergeron! (J2010J) pointed 
out, the strengths of some He I lines are especially gravity sensitive, and we make use of He I A3820, A4388, A4713 
and A4922, in that regard. We find the set of models: T e// /log g = 14000/8.1, 14400/8.4 and 14800/8.7, all give 
reasonable and roughly similar fits for most of the He I lines, shown plotted by the red line in Figure [2j It turns 



3 While G241-6 does have a parallax measurement I van Altena et al.lll995f ). it is not precise enough to improve on the spectroscopically 
assumed T e ///log g used bv lZuckerman et al. I |20ld ). 
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out that the over-prediction of He I A4713 mentioned bv lFriedrich et al.l (|1999l ) is resolved by our use of higher than 
average gravity, although there remains an unaccounted-for shift of that line's center (top panel Figure [2]). 

Another reason we prefer models with a higher than average gravity, is we noticed that at the value log g=8.0, for 
a given element, even though we could get a reasonably good fit to the equivalent widths with the same abundance, 
for almost all lines the model profiles are too deep and too narrow; this appeared to be the case for most elements. 
Having worked with matching similar models to similar HIRES observations of similar WDs (GD40 and G241-6), 
we presume that a mismatch in high-Z line broadening between the data and instrument-convolved model, is not a 
problem of accounting for the instrumental broadening, but rather is indication that there is some intrinsic broadening 
needed in the model, such as a higher gravity. We found that the models listed in Table [3] produced better overall 
fits to the high-Z line profiles than models of average gravity. Similarly, the fit to the Ha profile (Figure [3]) is better 
with higher gravity. 

Finally, with our choice models, the ionization balance for abundances of Mg and Fe naturally come out of the 
analysis of HS2253 at our desired consistency of within 0.1 dex between abundances from neutral and singly ionized 
states (Tabled]). Obtaining agreement in the ionization balance has proven challenging in other stars, including 
GD362, GD40 and PG1225 as described below. 



3.2. PG 1225-079 T e// /log g 

At temperatures < 12,000 K, the optical He I lines begin to disappear; He I A5876 (Figure [8]) is the strongest line 
observed. In this range of T e jj, only a very small fraction of the photospheric helium is involved in the formation of 
optical helium lines, since they arise from levels of high excitation energy, e.g. 20.9 eV above ground. Small changes 
in the atmosphere may have a large impact on the He line appearance, and we know that there is some unmodeled 
activity in this stars' upper atmosphere as is apparent from the emission cores of the Ca II H&K lines (Figure fl~5"j) 
and the narrow core of He I A5876 (Figures [8l &: 1 16[) . For the discussion of atmospheric parameters that follows, we 
ignore the narrow core of the He I A5876 feature described in Section [2j and only consider the broad component of 
He I A5876. 

The SEDs of our current models for PG1225, scaled to 2MASS photometry, are in good agreement with the 
IUE spectrum an d GALEX photometry (upper panel Figure [T7|) for T e yj near 10500 K, as used in previous work by 



Wolff et al.l (|2002l ) with log g = 8.0. However, looking at the ionization balance, a 10500/8.0 model results in a 0.2 dex 
difference between iron abundances derived from Fe I and Fe II. This is in the sense that Fe I lines are over-predicted 
and the discrepancy increases to 0.4 dex with a higher gravity model at 10500/8.3. Inversely, we find at T e ^y=10500 
K, the two states of Fe are brought into marginal agreement (within 0.03 dex uncertainties - Table [5]) by lowering the 
gravity to log 5=7.7. At this point in the T e ff/logg parameter space, the 10500/7.7 model produces a He I A5876 
line with a predicted EW that matches the measured value at 1.5 A, but the predicted profile is somewhat narrower 
and deeper than observed (red curve in Figure [8]). Also, the He I A4472 is not matched. The model helium lines can 
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be made shallower by increasing the gravity or decreasing the temperature, but always at the expense of degrading 
the ionization balance and the EW fit. As was the case for HS2253, there is a degeneracy of model solutions that 
fit the equivalent width of the helium line - for example 10500/7.7, 10800/8.0, and 11100/8.3 all produce a similar 
strength (and depth) He I 5876 feature. Accepting the imperfections in the helium line profile fitting, we favor these 
models as our best compromise since they: 1) are in reasonable agreement with the UV/optical/IR spectral energy 
distribution, 2) fit the He I A5876 equivalent width, and 3) support Fe I/Fe II abundance agreement to within 0.15 
dex. We use this set of three models for establishing pollutant-to-pollutant abundance ratios, which we show in 
Section^ are not very sensitive to the particular choice of model. 



4. ABUNDANCES 

Abundances are calculated for each line by linearly scaling the model input abundance by the ratio of EWs 
measured in the observed spectrum (Table [2]) to that of the model. In cases of blended features or otherwise noisy 
regions, a visual inspection of the profile fit aids the determinations. For each element we use the average of the set 
of calculated line abundances as the input for the next iteration of the model. It takes a few iterations until the 
average converges. Individual abundances are given equal weight in the final average element abundance. The EW 
measurement uncertainties are propagated through to an uncertainty in the final average abundance. We compare 
that to the mean deviation for the set of lines of an element, and assume the larger error. 

A ch allenge to abundan ce accuracy can come from discrepancies in the atomic line data and broadening pa- 



rameters (jVennes et al.ll201ll ). However, with multiple lines of each element observed, biases due to errors in the 
atomic data are minimized. We typically find a low dispersion in derived abundances for elements with multiple line 
detections (ranging from 0.02 - 0.15 dex) - generally comparable to, or smaller than, the propagated measurement 
uncertainties - suggesting that the errors due to line data (including broadening parameters) are not a big factor in 
the current results. 

Tables [5] and [5] list the pollutant abundances relative to helium, derived from the different models considered 
for each star. This explicitly shows that the choice of model temperature and gravity are a significant source of 
systematic uncertainty in absolute abundance determinations. The absolute (pollutant-to- helium) abundances vary 
by 0.2-0.4 dex with the modest range of model changes we apply. Also, estimates of the minimum accreted mass 
span a factor of ten, due to the fact that the calculated size of the convective zone has a strong dependence on 
the choice of model parameters as shown in Table El Fortunately for our goals, the relative (pollutant-to-pollutant) 
element abundances are rather homologous from one model to another with ratios varying only ~40% or less, a 
comparable contribution to the uncertainty as from the propagated measurement error. The upshot is that the 
relative composition of the polluting elements can be measured to a greater precision than that of the absolute mass, 
or rate, of material accreted. The relative element abundances with respect to the fiducial - magnesium - are given 
in Tabled) 
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Comparisons of abundances with prior work is possible for H, Mg, Si, Ca and Fe, although previous uncertainties 
are typically large (0.2-0.7 dex), and tend to propagate to a sizable range for the abundance rati os. For PG1225, th e 
HIRES abundances relative to helium are nearly all within the uncertainties of those quoted bv lWolff et al.l ([2002). 
Within the uncertainty limits, pollutant-to-pollutant ratios agree between previous and present analyses. From the 
HIRES data of PG1225, we obtain abundances for H, Mg, Ca, Sc, Ti, V, Cr, Mn, Fe, and Ni, and upper limits 



on Q , Na, Al, and Sr, with an upper limit for Si that is consi stent with its UV detection and analysis (jWolff et al 



2002]); an upper limit for carbon is taken from the UV work (IWolff et all 120021). For H S2253. the nominal HIRES 



abundances with respect to helium are all lower than those quoted by iFriedrich et al.l (|1999l ) by at least 0.6 dex, 
most likely due to the differences in T e ff/logg model inputs, particularly the much larger gravity used in the pre sent 
work. Nonetheless, the pollutant-to-pollutant abundances derived from HIRES data and IFriedrich et al.l (|1999l ) are 
mostly in agreement, with the exception of the Fe/Mg ratio. Previously, HS2253 appeared iron-rich (by number, 
Fe/Mg = 2, with 1.2 to 3.3 covering the range of uncertainty); the present HIRES analysis yields a ratio that is 
nominally half as large (Fe/Mg = 0.89 ± 0.18), and similar to Solar and Bulk Earth values. We obtain abundances 
in HS2253 from the HIRES data for H, O, Mg, Si, Ca, Ti, Cr, Mn, Fe, and a marginal detection of Ni; we derive 



uppe r limits on Na, Al, Sc, V, and Sr, while an upper limit for carbon is taken from the UV work (IFriedrich et al 
1999h . 



5. ACCRETION DIFFUSION MODELS 

As mentioned in Section [1] a helium-dominated polluted WD system has an evolving interplay of accretion and 
diffusion, since these WDs have relatively long settling times. Calculated settling times specific to our p articular 



model s of HS2253 and PG1225 for the various elements studied are given in Tables 0] and [5j As described in lKoester 



(2009), the accretion-diffusion situation is roughly categorized by 3 phases: 1) an early phase in which accretion 
began "recently", that is, before differential settling has had time to significantly alter the abundance ratios - the 
observed atmospheric abundances directly represent the accreted ones, 2) a steady state with ongoing accretion, 
where the observed abundances are related to the accreted ones by ratios of settling times, and 3) a dissipation 
phase after accretion has slowed, paused or ended, and elements with shorter settling times disappear faster than 
those with longer settling times. We use these benchmarks to characterize our observations. 

If accretion has been going on longer than a couple of settling times, then a steady state interpretation may 
be appropriate. This assumes a semi-continuous flow of material lasting for ~1 Myr or more in the systems under 
study here. To calculate the steady state abundances, for each element we take the average over different models of 
the ratios of settling times, and use the averag e ratios to app ly a differential settling "correction" to the abundance 



ratios of Table [6] as prescribed by Equation 7 of lKoesterl (|20091 ) . It is worth noting that analogous to the abundances, 



while the settling times themselves vary significantly from one model to another, the pollutant-to-pollutant ratios of 
settling times vary very little (RMS of most < 3%, all < 8%) with the model variations for a given source. For each 
element, the RMS's of diffusion data among models are propagated in quadrature into the steady state abundance 
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uncertainties. 

The principle result of applying a steady state "correction" is an increase in the abundance of iron and other 
heavy elements relative to the lighter elements in the inferred parent body. This effect is apparent in Figure [18] where 
the early phase (= connective zone = CVZ) abundances are plotted as filled markers, while the the open markers 
represent parent body abundances for the system in a steady state. For the set of elements plotted in Figure [TH] Mg 
and Si are the lightest and have the longest settling times. Since Mg is the fiducial, applying a settling "correction" 
moves all the points of the plot up, relative to Mg, as seen in going from early phase to steady state values. The 
difference between an early phase and a steady state interpretation is most pronounced in the O/Fe ratio, which can 
be seen by comparison of the bulk composition bar graphs, Figures \19\ and 1201 



6. RESULTS AND DISCUSSION 



From Table [6] we find the accreted abundances are not in Solar proportions for the volatiles H, C, and O 
(HS2253), verifying that these systems are neither accreting from the interstellar medium nor gas giant planets. 
Rather for example, the distinctively low carbon abundances in both HS2253 and GD40 are such that C/Mg is more 
than an order of magnitude below the Solar abundance and therefore more similar to Earth (Jura 2006). It is helpful 
to compare all measurements to a set of values such as the bulk Earth. In the following we will discuss the abundance 
ratios of Table El mainly through Figure [TSl in which the observed abundance ratios are normalized to Mg and the 
bulk Earth, and Figures [T9l and 1201 in which we compare the major element compositions of inferred parent bodies 
with the bulk Earth. 



HS2253 and PG1225 eac h have the high- Z mass of a large as teroid ~10 g cu rrently in their convective zones 
(Table EJ - similar to GD362 (jZuckerman et al.ll2007h and GD40 ((Klein et alJboioh . HS2253 has one of the highest 



time-averaged mass accretion rates of all polluted white dwarfs, being comparable to GD362 (jjura et al.l 12009 ; 



Koester 



2009|) at 



10 10 g 8 _1 , 



with SDSS0738 apparently bein g the cu r rent champion at ~10 g s based on 



Dufour et al.l (|2010l ) and our settling time calculations parallel to iKoesterl (|2009l ). The time-averaged mass accretion 
rate of rocky material in PG1225 is a few xlO 8 g s _1 , about a factor of 5-10 less than GD40. 



6.1. HS 2253+8023 Composition 



If a WD accretes a rocky planet ary body, then the relative element abundances will have the distinctive propor- 



tions that arise from mineral oxides (jJura et al.ll2009l ; iKlein et al.ll2010l ). made up of MgO, AI2O3, Si02, CaO, TiO, 
Ci"203, etc. with iron coming from either FeO, Fe203, or Fe metal if a differentiated core formed. Another possible 
vehicle for oxygen in a terrestrial-type planetary body is water ice, which depends on how much hydrogen is present. 
For HS2253, the HIRES measurement of oxygen plus all major elements, some minor ones, and hydrogen, allows us 



to evaluate the balance of oxides and water in this system according to Equations 2 and 3 of 



Klein et al. 



ipoich . 



Looking at the oxide balance in an early phase interpretation for HS2253, we find that the number ratio of 
oxygen that can be accounted for by oxides with the other elements (using half the iron in FeO, half in Fe2C>3, 
plus aluminum at an Earth-like ratio of Al/Mg = 0.086) is O/Mg = 3.8 ± 0.6. Since the observed amount of 
oxygen is O/Mg = 5.4 ± 1.1, at the limits of the uncertainties the oxygen can be carried in rocky minerals if iron is 
predominantly in oxide form; water is not needed to balance the oxides. Generally though, due to the ambiguity of 
the state of iron and since iron is a major element, there is a sensitive trade-off between metallic iron and water in 
balancing the oxides. In HS2253 the early phase abundances can be attributed solely to mineral oxides, but could 
involve up to ~15% water by mass (limited by the observed hydrogen abundance). 

In the steady state for HS2253, the number ratio of oxygen expected to be associated with the observed major 
elements, Mg, Al, Si, Ca, and Fe, as minerals (half the iron in FeO, half in Fe2 03, plus an Earth-like ratio for Al) is 
O/Mg = 4.7 ± 0.8, and the ratio of of observed oxygen coming from the parent body is O/Mg = 4.5 ± 0.9. In the 
steady state the ratios agree very well with a rocky mineral oxide composition. Again, since some Fe may have been 
metallic and there is a range of uncertainty overlap in the O/Mg ratios, there is still sufficient oxygen to allow for 
some water; the amount of possible water in the parent body under a steady state interpretation is by mass <7%, 
limited by the observed hydrogen abundance. While HS2253 does not display an infrared excess, it can be seen from 
Figures [19] and [20] that this heavily polluted WD has atmospheric abundance ratios that are similar to GD40 - an 
analogous WD that does display an infrared excess. 

Could the system be in a dissipating phase? Not likely, since the oxides and water are well balanced in both 
the early phase and steady state interpretations; a dissipative interpretation will eventually run into difficulty with 
oxides. From the lower panel of Figure [18] in an early phase interpretation, the abundances of Ti, Mn, Cr, Fe, and 
Ni, are in excellent agreement with bulk Earth values (the Ca overabundance and Si underabundance are similar to 
GD40). As accretion-diffusion proceeds, the ratios move up to the steady state values. If dissipation takes over, the 
inferred abundance ratios would be driven even further up compared to Mg and Si, while O/Mg decreases, and would 
imply a parent body that is somehow depleted of the major elements Mg, Si, and O, compared to high abundance 
ratios of Ca, Ti, Mn, Cr, Fe and Ni. We have no physical explanation for how that may occur. Thus, we prefer 
the simpler explanation that HS2253 most likely accreted material with abundances predominantly similar to bulk 
Earth and is somewhere at or between an early phase and a steady state of accretion and diffusion. 



6.2. PG 1225-079 Composition 



In addition to significant high-Z element pollutio n, PG1225 has a considerable amount of hydrogen in its 
convective zone. Where did it come from? For GD 362, lJura et al.l (120091 ) showed that the presence of circumstellar 
dust, high-Z mate rial and copious h ydrogen can simultaneously be explained by the accretion of a large ice-rich 
planet. Figure 1 of I Jura et al.l (J2009J) illustrates that the atmospheric pollution by hydrogen in GD362 and GD16 is 
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anomalously high compared to other known helium-dominated externally polluted WDs 4 of their evolutionary epoch, 
10,500 K < T e jf < 11,500 K, i.e. cooling age 0.4 - 1 Gyr. At a similar cooling age, our set of possible models for 
PG12 25 indicate a hydrogen mass which can range from6xl0 22 to 8xl0 23 grams. Referring to Figure 1 of lJura et al.l 



(2009) this is similar, perhaps on the high-end, but not extraordinary, compared to the average amou nt of accreted 
hydrogen in the set of helium-dominated WDs with detected hydrogen at this age. lJura &: Xul (120101 ) suggest that 
the origin of the hydrogen in these more ordinary WDs may be from the accretion of ice-rich planet (esimal)s which 
survive post-main sequence evolution. 

Unfortunately, we cannot evaluate a water and oxide balance for PG1225 since the Ol A7775 lines are not visible 
in its HIRES spectrum. If the detected elements came from a rocky parent body, then the amount of oxygen (by 
number) associated with mineral oxides is O/Mg ~ 5, which is far below the sensitivity of our observations that yield 
an upper limit of O/Mg < 60 (Tabled]). With just a small portion of the allowed O/Mg number ratio allocated to 
rocky oxides, the remaining number ratio that could be associated with water is O/Mg < 55, but that only accounts 
for less than 10% of the observed hydrogen. In order to explain all of the observed hydrogen as being derived from icy 
parent bodies, either an ongoing steady state of atmospheric settling from a somewhat large parent body (perhaps 
a Ceres analog) or previous ice-rich accretion event (s), would be implied. 

PG1225's spectrum is very rich with high-Z lines and the numerous newly detected elements display an interest- 
ing abundance pattern. In the top panel of Figure [181 Mg, Cr, Mn, Fe and Ni, have abundance ratios similar to bulk 
Earth values, but we find that the refractory elements - Ca, Sc, Ti and V - are all a factor of ~2-3 higher abundance 



than in the Earth . High Ca and T i abundance ratios have previously been measured in GD362 (jZuckerman et al 



20071 ) and GD40 (IKlein et al.ll2010l ). As shown in the bottom panel of Figure [THl HS2253 also has a somewhat en- 
hanced Ca/Mg ratio. The enhancement of refractory element abundances with condensation temperatures >1400K 
is significant and uniform in PG1225, and the upper limits derived for Al and Sr are consistent with the trend. 
One possible explanation for this observation is that the material accreted onto PG 1225-079 derived f rom a plan 



etary body that formed in a higher temperature environment than did the Earth. The simulations of iBond et al 



()2010l ) suggest that refractory-rich exoplanets should be a common outcome of planet formation in nebular regions 
sufficiently close to the host star. 



PG1225 has ju st a subtle infrared excess appearing only at 7.9/im (jFarihi et al.l l2010bl ; iKilic et al.l l2008al ). 
Farihi et al.l (|2010bl ) raise the possibility that it could be due to a relatively narrow dust ring with an inner hole cleared 
by accretion, suggesting that the system is in a dissipative phase. Referring to Figure [TBI since all the refractory 
elements have shorter settling times than Mg, any "correction" for settling causes an even greater overabundance of 
refractory species inferred in the parent body. In other words, contrary to the expectation from a dissipative phase, 
the heavy elements do not appear sunk. However, the settling times in this WD are quite long ~Myr, and disk 



4 See also iTremblav et al.l ( 201 if ) for 2 more recently uncovered candidates of this unusual type of WD - helium-dominated with large 
quantities of atmospheric hydrogen. 
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lifetimes may be shorter; a large range has been suggested ~10 5 -10 9 yr ( Jura et alj 12009 . and references therein). 
Thus it is conceivable that we could be viewing the system at a point toward the end of an accretion episode but 
before dissipative differential settling has a significant effect on the atmospheric ratios. Whatever the case regarding 
circumstellar dust, the atmospheric element ratios imply that PG1225 most likely accreted material with abundances 
predominantly similar to bulk Earth, but enhanced in refractory species, and the system is somewhere in or near an 
early phase or steady state of accretion and diffusion. 



6.3. Helium-Rich Highly-Polluted White Dwarfs 



Finally we discuss the set of highly polluted helium-dominated white dwarfs that have been comprehensively 
studied. To date, 5 helium-dominated WDs have >8 elements reported (GD362, GD40, G241-6, HS2253+8023, 
PG1225-079) plus 2 more (SDSS0738, GD61) have measurements of all major (O, Mg, Si, and Fe) terrestrial planetary 
constituents. 

Though our focus in this paper is on compositions in heavily polluted, helium-dominate d WDs, we recognize the 



importance of the study of heavily polluted, hydrogen- dominated WDs, such as NLTT43806 (IZuckerman et al 



and GALEX J19 3156.8 + 011745 (jVennes et al.ll2010h . With the detection of circumstellar dust (jDebes et al 



2011) 



2011 



Melis et allEoil] ) and all the major element species (O, Mg, Si, and Fe), GALEX J193156.8 + 011745 (GALEX1931) 
is almost certainly accreting a planetary body; with settling times on the order of weeks, accretion must be happening 
now. However we note there is curre ntly a troubling inconsistency between the modeled results of different groups 



(Vennes et al 



2011 



Melis et al.ll201ll ). which is due in part to a large, unexplained difference for the assumed effective 



temperature (AT e jj > 2,000 K) of GALEX1931, and in part to the treatment of atmospheric settling, which for 
warm (T e ff ~ 20,000 K) hydrogen-dominated WDs, is complicated by the lack of a convective envelope. 

Figures [19] and [20] represent the compositions by mass of major elements in the helium-dominated polluted WDs 
in which oxygen has been detected, or a meaningful upper limit set (GD362). Lacking an oxygen detection, PG1225 
is not on the bar graphs. There are two different bar graphs showing two different interpretations of the parent 
body abundances, in the early phase and in the steady state. The graphs are related by ratios of settling times as 
described in Section [5] We are interested in comparing the bulk accreted compositions, and since hydrogen never 
settles out of the convective zone, we cannot be certain how much of it is associated with the current heavy element 
pollution. We do not include hydrogen in Figures [19] and 1201 which has a negligible effect on the graph appearance 
for bulk Earth, HS2253, GD40, SDSS0738 and G241-6, being less than 1% of the accreted mass in those sources. 
On the other hand, in GD362 hydrogen is over 99% of the accreted mass and in GD61 it is over 75%. As discussed 
in the preceding section, scenarios have been proposed to explain the presence of hydrogen from the accretion of 
ice-rich planetary bodies. Now we consider the compositional possibilities from the perspective of oxygen. 



In each of the WD systems with detected oxygen, there is always enough O to account for the delivery of the 
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major 5 elements Mg, Si, Ca, and Fe, in the form of a terrestrial-type planetary body. The upper limit of oxygen in 
GD362 is just enough to have come from Mg, Si, Ca and Fe in oxides, with the possibility of metallic iron. Since there 
is very little "spare" oxygen in GD362, in or der to acco u nt for the hydrogen, high-Z elements, infrared excess, and 
lack of x-ray luminosity in a unified picture, lJura et al.l (|2009l ) calculate a scenario in which settling and accretion 
have been ongoing for a long time, implying the consumption of a large ice-rich planetary body at least the size of 
Callisto. 

On the right side of the bar graphs, the abundance ratios of GD61 imply the accretion of more oxygen than 
can be carried in rock. An O "excess" can be the result of dissipative gravitational settling, which would make 
the abundances seem oxygen-rich and i ron-poor. Anothe r way to get an O excess is by the accretion of an ice-rich 
parent body, as considered in detail by I Jura Xul (120101 ). The especially high O/Fe and Mg/Fe abundance ratios 
in the atmosphere of GD61 look suspiciously similar to what is expected in a system where accretion has stopped or 
paused - the heavier elements sink faste r than the lighter ones. However, along with measurements of all the major 
element abundances, iFarihi et al.1 (|201ll ) report an infrared excess for GD61. Sinc e the W D has plenty of hydrogen 
in its atmosphere, an O excess, and detected circumstellar material, IFarihi et al.l (|201ll ) raise again the possibility 
of the accretion of an ice-rich planetary body. Due to abundance measurement uncertainties and the uncertainty of 
the accretion-diffusion phase, the results are yet inconclusive regarding the composition of the parent body accreted 
onto GD61. 

G241-6 also appears to be a case with an O excess, even at the limits of the abundance uncertainties. The upper 
limit on hydrogen in the atmosphere amounts to less than 0.6% of the mass of detected contaminants, and therefore 
cannot sufficiently account for the excess oxygen by an origin in ice; the parent body must have been relatively dry. 
The abundance pattern of G241-6 is similar to GD40 but with Mg, Si, Ca and Fe appearing "sunk" (compare the 
GD40 early phase bar of Figure [19] with G241-6's steady state bar of Figure [20]) . In other words, starting from GD40 
abundances, if we could turn accretion off, then eventually gravitational settling will evolve the abundances into 
ratios similar to what we observe for G241-6. Since we now know that G241-6 does not have an infrared excess (Xu 
& Jura, Spitzer data, in preparation), it is likely that G241-6 is in a dissipative phase, having originally accreted an 
ice-free planetary body with composition similar to GD40 and the bulk Earth. 

The abundances in the remaining three: HS2253, GD40, SDSS0738, are all consistent with the major elements 
being derived from rocky mineral oxides, and a composition similar to bulk Earth being composed more than 85% 
by mass in O, Mg, Si, and Fe, for either a steady state and/or early phase interpretation. Overall their parent body 
compositions are depleted of volatiles; in addition to the dearth of carbon in GD40 and HS2253, all of them have 
distinctively low abundances of hydrogen (compared to high-Z elements) and little, if any, excess oxygen - again 
indicating a dearth of ice in the accreted material. 



5 Aluminum may also be expected in comparable mass as Ca, but it is more difficult to detect. Within uncertainties, the undetected 
presence of Al at an Earth-like abundance will not substantially change the oxide balance. 



- 15 - 



While GD362 and GD61 are candidates for exoplanetary systems having accreted a substantial amount of water, 
the other four stars with oxygen measurements must have accreted relatively dry planetary bodies (<10% H2O by 
mass). Furthermore, in all four of those cases, the pollution is most plausibly dominated by planet (esimal)s that 
formed with little or no ice, for reasons as follows. GD40 and SDSS0738 both have an infrared excess, and as 
discussed by Jura (2008), the existence of a dust disk is best understood as resulting from the disruption of a single 
parent body, rather than multiple parent bodies - which, arriving on trajectories of varied inclinations, tend to result 
in gaseous disks due to dust grain destruction by sputtering in high- velocity impacts. With a minimum accreted 
mass of ~10 22 g, and a rock density o f 3 g cm" 3 , the ra dius of a single polluting parent body is greater than ~100 



km. Thermal erosion calculations by lJura fc Xul (|2010l ) demonstrate that asteroids of radii > 100 km would have 
retained at least 50% of their internal ice, if present, through red giant evolution. Thus the parent bodies polluting 
GD40 and SDSS0738 most likely formed without much ice. 

Could HS2253 and G241-6 be polluted by many small (radius < 20-30 km), initially ice-rich asteroids whose 
internal water was lost during the host star's Asymptotic Giant Branch (AGB) phase? Consider a model for HS2253 
involving 2 5 km a steroids that were 25% water by mass before the star was on the AGB. According to Figure 4 of 



Jura k, Xul (|2010l ). unless they originally orbited at less than 5 AU, 25 km asteroids would retain at least 10% of 
their water during AGB evolution of a 3 M star. At a density of 2.1 g cm -3 (analogous to Ceres), the pre- AGB 
asteroid might have had 10 20 g of high-Z material, and 4 x 10 18 g of hydrogen (H) from ice. After the AGB, the 
asteroid has 10 20 g of high-Z material, and 4 x 10 g of H. With the middle model from Table (H the convective 
zone of HS2253 has ~2 x 10 22 g of high-Z matter. Therefore, it needs to have accreted about 200 asteroids during 
the most recent 10 5 years, or about 1 asteroid every 500 years. During this time, the white dwarf would acquire 8 
x 10 19 g of H, about a quarter of the total H (~3 x 10 20 g with the middle model from Table |4|), which the WD 
has accumulated over its entire cooling age of approximately 400 Myr. The relatively small amount of H in the 
convective zone implies that no more than ~600 similar asteroids could have arrived in the previous 400 Myr. Why 
would the asteroids arrive at a rate of 0.002 yr _1 recently and 2 x 10~ 6 yr -1 previously? The result of this model 
exercise is analogous, but more extreme, for G241-6. Another difficulty with the above hypothesis is that, even if 
there is a burst of impacts, one imagines that there would be a size distribution of impactors and, at least in the usual 
kind of model, most of the mass would be carried in the larger objects. All that is required to produce the observed 
hydrogen is the accretion of single ice-rich asteroid of radius 100 km. Unless the size distribution of asteroids is very 
steep, there could easily be one big asteroid for every 200 small ones. Therefore, to explain the pollutions of HS2253 
and G241-6 by swarms of initially ice-rich small asteroids, would require both a spike in time of number frequency 
and and a spike in the size distribution of the asteroids. This double constraint seems implausible. It is most likely 
that the planet (esimal)s which pollute these four now comprehensively studied WDs - HS2253, GD40, SDSS0738 
and G241-6 - were formed with little ice, and we are measuring the distilled elemental compositions of bodies which 
condensed out of protoplanetary disks at radii interior to any snow line. 
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7. CONCLUSIONS 



We have obtained and analyzed high-resolution optical spectra of HS2253+8023 and PG1225-079 with HIRES 
at Keck Observatory. This work adds two more systems to the small but fast-growing sample of polluted white 
dwarfs with detailed and comprehensive analyses. Four out of six objects with oxygen measurements most likely 
accreted planet (esimal)s of mineral oxide constitution, which were nearly ice-free, consistent with formation interior 
to a snow line. It is clear that extrasolar planetary systems produce rocky bodies that are compositionally similar to 
terrestrial planets in our own solar system; Earth-like planet (esimal)s apparently do form elsewhere in the galaxy. 
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Table 1: Keck/HIRES Observations 



7"P\ 


TTrp I "\ i 

U 1 Date 


Approx. Seeing 
(visual) 


Collimator 


A range 
(A) 


Resolution 

A 

AA 


Exposure 
(seconds) 


C OA 1 R 


ZUU / IN ov IX) 


rY'a 
U. o 


-Drue 


q 1 on ^ot^n 


4UUUU 






2008 Aug 06 


l'/3 


Blue 


3140 - 5950 


40000 


2 x 1500 




2008 Nov 14 


l'/l 


Red 


4500 - 9000 


40000 


1800+2000 


HS 2253+8023 


2008 Aug 06 


l'/3 


Blue 


3140 - 5950 


40000 


1800+2400 




2008 Nov 14 


l'/l 


Red 


4500 - 9000 


40000 


2 x 2000 




2008 Nov 15 


l'/4 


Red 


4500 - 9000 


28000 


2 x 1800 


PG 1225-079 


2008 Feb 13 


0"8 


Blue 


3120 - 5950 


40000 


2 x 1600 




2008 Feb 14 


not avail 


Blue 


3120 - 5950 


40000 


2600 




2008 Feb 26 


0"8 


Red 


4600 - 9100 


40000 


4000 



On 2008 Feb 14 the sky was partly cloudy with an estimated visual extinction of ~0.5 mag. The rest of the 
observations occurred during mostly clear skies. Gaps in blue coverage with our setups occur between 4010-4100 A, 
and 5990-6100 A. There are many small gaps in red coverage. 



Table 2:: Absorption Line Measurements of HIRES Spectra 



Ion 


A a 


X 


log gf 




Equivalent Width (mA) 






(A) 


(eV) 




PG 1225-079 


HS 2253-8023 


G241-6 c 


H I 


6562.79 


10.20 


0.71 


10,500 (2,000) 


850 (120) 


< 250 


H I 


4861.32 


10.20 


-0.02 


3,500 (1,500) 


120 (30) 




O I 


7771.94 


9.15 


0.37 


< 250 


350 (90) 


250 (50) 


O I 


7774.16 


9.15 


0.22 




190 (50) 


170 (50) 


O I 


7775.39 


9.15 


0.00 




220 (110) 


90 (40) 


O I 


8446.35 


9.52 


0.24 




380 (190) 


210 (50) 


Nal 


5889.95 





0.12 


< 50 


< 120 




Mg I 


3829.35 


2.71 


-0.23 


in 3800 interval 


50 (6) 


15 (5) 


Mg I 


3832.30 


2.71 


-0.36, 0.12 


in 3800 interval 


220 (40) 


75 (14) 


Mg I 


3838.29 


2.72 


0.39, -0.35 


in 3800 interval 


470 (40) 


106 (14) 


Mg I 


5172.68 


2.71 


-0.45 


30 (15) 


70 (40) 


31 (15) 


Mg I 


5183.60 


2.72 


-0.18 


90 (45) 


170 (40) 


52 (11) 


Mg II 


4481 


8.86 


0.74,-0.56,0.59 


90 (20) 


1040 (110) 


470 (50) 


Mg II 


7877.05 


10.00 


0.39 




270 (110) 


110 (25) 
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Table 2 — Continued 
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Table 2 — Continued 
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Table 2 — Continued 
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Table 2 — Continued 
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64 (17) 




Fe II 


3258.77 


3.89 


-1.15 


10 (3) 


48 (13) 




Fe II 


3259.05 


3.90 


-0.97 


16 (5) 


68 (24) 




Fe II 


3277.35 


0.99 


-2.39 


20 (6) 


66 (14) 




Fe II 


3281.29 


1.04 


-2.74 


13 (4) 


30 (8) 




Fe II 
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3.81 
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3323.06 
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Table 2 — Continued 
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losi a f 


Equivalent Width (mA) 
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rlo ZZDo-oUZo 
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2.58 


-2.51 




21 (9) 




Fe II 


4233.17 
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Fe II 
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Fe II 
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2.84 


-2.25 
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Fe II 


4549.47 


2.83 
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Fe II 
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Fe II 


4923 93 


2.89 
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DO (14) 


77 (15) 




Fe II 
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1 /in / 1 O \ 

142 (13) 
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no / A CT \ 
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21U (20) 
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Fe II 
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on / rr\ 
oU (0 ) 




Fe II 


5316.61 


3.15 


-1.87 




70 (29) 




Ni I 


3414.76 


0.03 


-0.014 


30 (12) 






Ni I 


3524.53 


0.03 


0.008 


40 (16) 






Ni I 


3619.39 


0.42 


0.04 


in Fe I 3618.77 






Ni II 


3513.99 


2.87 


-1.46 




18 : 




Sr II 


4077.71 





0.17 


< 60 


< 90 




interval 


3226-3231 






450 (5Q)b-FeII+TiU 


see Fell, Till 


see Fell 


interval 


3800-3850 






2,100 (400) 6 - M ^+^e/ 


see Mgl, Fel 


see Mgl 



a Laboratory wavelength in air. Wavelength, x (excitation energy of lower state) and log gf data from the Vienna 
Atomic Line Database. 

b Blended line - contributing element noted in italic with line data given in the table (Table [2]) . The equivalent width 
is from the total summed fl ux of the blend. 

c G241-6 abundances from IZuckerman et al.l (|2010l ) are associated with the EWs of this table. Spectra of G241-6 
are shown in Figure [13] of this paper and Fig ure 4 of IZuckerman et al.1 (120101 ) . 

d Combined Ca H&K EW measurement from Koester et al. ( 2005bl ) is used since the HIRES data do not fully cover 
the red wing of this large feature. 
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Table 3: White Dwarf Models, Derived Masses and Accretion Rates 





T eff 


log g 


M* 


log M^/M* 


thigh-'/, 


< M acc > 




(K) 


(cgs) 


(M ) 




(10 22 g) 


(g s- 1 ) 


PG1225 model 1 


10500 


7.7 


0.42 


-4.30 


5.6 


3 x 10 8 


PG1225 model 2 


10800 


8.0 


0.58 


-5.07 


2.1 


5 x 10 8 


PG1225 model 3 


11100 


8.3 


0.77 


-5.84 


0.7 


9 x 10 8 


HS2253 model 1 


14000 


8.1 


0.65 


-5.77 


6.0 


5 x 10 9 


HS2253 model 2 


14400 


8.4 


0.84 


-6.48 


2.0 


9 x 10 9 


HS2253 model 3 


14800 


8.7 


1.03 


-7.21 


0.6 


1.4 x 10 10 



Convective zone (CVZ) helium masses, contaminant (high-Z) masses, and time-averaged mass accretion rates (ex- 
cluding hydrogen) for WD models of varying effective temperature (T e ff) and gravity (g). The (steady state) 
accretion rates are the sum of individual element rates calculated from the mass of the element in the CVZ divided 
by its settling time, with abundances and settling times from Tables 2] and [5j For PG1225-079 the Si upper limit is 
used, and an abundance of oxygen is assu med as would deriv e from the dominant mineral oxides: MgO, Si02, FeO, 
CaO. Including carbon at its upper limit (jWolff et al.ll2002l ) increases the high-Z masses and accretion rates by ~ 
30% for PG1225-079, and is negligible for HS2253+8023. Including aluminum at its upper limit increases the high-Z 
masses and accretion rates by ~5% for both stars. 
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Table 4: HS2253+8023 Atmospheric Abundances by Number and Settling Times for 3 Different T e yj/log g Models 



z 


14000/8.1 
log(Z/He) 


14400/8.4 
log(Z/He) 


14800/8.7 
log(Z/He) 


14000/8.1 
log(r/yrs) 


14400/8.4 
log(r/yrs) 


14800/8.7 

log(r/yrs) 


H 


-5.70 (0.06) 


-5.62 (0.06) 


-5.55 (0.06) 











-5.48 (0.07) 


-5.37 (0.07) 


-5.28 (0.07) 


5.71 


5.01 


4.31 


Mg 


-6.23 (0.06) 


-6.10 (0.04) 


-6.02 (0.03) 


5.63 


4.93 


4.23 


from Mg I 


-6.27 (0.05) 


-6.12 (0.06) 


-6.03 (0.04) 








from Mg II 


-6.17 (0.05) 


-6.07 (0.05) 


-6.01 (0.05) 








Si 


-6.40 (0.03) 


-6.27 (0.03) 


-6.16 (0.04) 


5.63 


4.93 


4.21 


Ca 


-7.10 (0.03) 


-6.99 (0.03) 


-6.92 (0.05) 


5.50 


4.78 


4.06 


Ti 


-8.92 (0.02) 


-8.74 (0.02) 


-8.55 (0.02) 


5.43 


4.73 


4.03 


Cr 


-8.17 (0.04) 


-8.01 (0.03) 


-7.83 (0.03) 


5.44 


4.74 


4.01 


Mn 


-8.63 (0.09) 


-8.42 (0.10) 


-8.24 (0.09) 


5.43 


4.71 


3.99 


Fe 


-6.34 (0.04) 


-6.17 (0.03) 


-6.00 (0.02) 


5.42 


4.71 


3.98 


from Fe I 


-6.29 (0.05) 


-6.15 (0.04) 


-5.99 (0.05) 








from Fe II 


-6.36 (0.03) 


-6.18 (0.03) 


-6.01 (0.03) 








Ni 


-7.51 : 


-7.31 : 


-7.14 : 


5.40 


4.68 


3.96 


C 








5.74 


5.05 


4.35 


Na 




< -6.8 




5.63 


4.92 


4.21 


Al 




< -6.7 




5.61 


4.92 


4.22 


Sc 




< -9.8 




5.45 


4.74 


4.03 


V 




< -8.6 




5.42 


4.73 


4.02 


Sr 




< -9.6 




5.24 


4.53 


3.81 



Element abundances and settling times derived from the equivalent width measurements of Table [2] and different 
models, as described in Sections El HI and [U Uncertainties in parentheses are either the standard deviation of the 
mean for the set of individual line abundances, or the propagated measurement error, whichever is greater. When 
a significant difference is measured between ionization states of an element (only true for iron), this additional 
uncertainty is added in. While the absolute abundances and diffusion times vary considerably between models, the 
relative abundances and relative diffusion times among the listed elements are far less sensitive to changes in the 
model. 
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Table 5: PG1225-079 Atmospheric Abundances by Number and Settling Times for 3 Different T e jj/log g Models 



z 


f0500/7.7 
log(Z/He) 


10800/8.0 
log(Z/He) 


11100/8.3 
log(Z/He) 


10500/7.7 
log(r/yrs) 


10800/8.0 
log(r/yrs) 


11100/8.3 

log(r/yrs) 


H 


-4.13 (0.10) 


-4.05 (0.10) 


-3.99 (0.10) 








Mg 


-7.44 (0.05) 


-7.27 (0.05) 


-7.11 (0.05) 


6.85 


6.14 


5.41 


from Mg I 


-7.44 (0.06) 


-7.27 (0.06) 


-7.11 (0.06) 








from Mg II 


-7.43 (0.10) 


-7.30 (0.10) 


-7.11 (0.10) 








Ca 


-8.25 (0.03) 


-8.06 (0.03) 


-7.90 (0.04) 


6.83 


6.07 


5.32 


from Ca I 


-8.30 (0.15) 


-8.18 (0.15) 


-8.04 (0.15) 








from Ca II 


-8.24 (0.03) 


-8.04 (0.03) 


-7.89 (0.04) 








Sc 


-11.57 (0.07) 


-11.29 (0.07) 


-11.05 (0.07) 


6.78 


6.03 


5.28 


Ti 


-9.68 (0.02) 


-9.45 (0.02) 


-9.23 (0.02) 


6.76 


6.01 


5.26 


V 


-10.68 (0.10) 


-10.41 (0.10) 


-10.18 (0.10) 


6.73 


5.98 


5.23 


Cr 


-9.47 (0.06) 


-9.27 (0.06) 


-9.05 (0.06) 


6.72 


5.97 


5.24 


Mn 


-10.03 (0.14) 


-9.79 (0.14) 


-9.52 (0.14) 


6.70 


5.96 


5.24 


Fe 


-7.64 (0.03) 


-7.42 (0.07) 


-7.21 (0.08) 


6.69 


5.97 


5.25 


from Fe I 


-7.67 (0.03) 


-7.50 (0.03) 


-7.30 (0.03) 








from Fe II 


-7.62 (0.03) 


-7.36 (0.03) 


-7.15 (0.03) 








Ni 


-8.88 (0.14) 


-8.76 (0.14) 


-8.64 (0.14) 


6.70 


5.99 


5.24 


C 








6.89 


6.18 


5.47 





< -5.69 


< -5.54 


< -5.44 


6.86 


6.15 


5.43 


Na 


< -8.45 


< -8.26 


< -8.02 


6.83 


6.12 


5.41 


Al 


< -8.04 


< -7.84 


< -7.68 


6.84 


6.11 


5.38 


Si 


< -7.49 


< -7.27 


< -7.07 


6.85 


6.11 


5.39 


Sr 


< -11.90 


< -11.65 


< -11.38 


6.56 


5.85 


5.10 



Similar to Table H but for PG1225-079. 
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Table 6: Abundance Ratios By Number 





PG 1225-079 
Convective Zone 


HS 2253±8023 
Convective Zone 


Solar 
CLodders 20031 


H/Mg 


1700 (500) 


3.1 (0.6) 


2.8 x 10 4 


C/Mg 


< 5 a 


< 0.01 b 


6.9 


O/Mg 


< 60 


5.4 (1.1) 


14 


Na/Mg 


< 0.12 


< 0.22 


0.056 


Al/Mg 


< 0.31 


< 0.30 


0.081 


Si/Mg 


< 1.1 


0.68 (0.10) 


0.98 


Ca/Mg 


0.16 (0.02) 


0.13 (0.02) 


0.062 


Sc/Mg 


9.4 ±2.7 x 10~ 5 


< 2.2 x 10~ 4 


3.3 x 10~ 5 


Ti/Mg 


6.7 ± 1.2 x 10~ 3 


2.4 ±0.6 x 10~ 3 


2.3 x 10~ 3 


V/Mg 


7.1 ±2.2 x 10" 4 


< 3.7 x 10~ 3 


2.8 x 10~ 4 


Cr/Mg 


0.010 (0.002) 


0.013 (0.003) 


0.013 


Mn/Mg 


3.1 ± 1.2 x 10~ 3 


5.1 ± 1.7 x 10~ 3 


8.9 x 10~ 3 


Fe/Mg 


0.71 (0.12) 


0.89 (0.18) 


0.83 


Ni/Mg 


0.033 (0.011) 


0.06 : 


0.047 


Sr/Mg 


< 4.5xl0" 5 


< 3.1xl0~ 4 


2.3 x 10~ 5 



a C/Mg from Wolff et al. (2002). 
b C/Mg from Friedrich et al. (1999). 

Measured atmospheric abundance ratios calculated as the average of 3 models with ratios from Tables U] and [5j 
Uncertainties in parentheses are a combination of the propagated abundance ratio measurement uncertainties, and 
the root-mean-square of ratio values from different models from Tables H] and El added in quadrature. 
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Fig. 1. — Combined, echelle-order-merged, shortest wavelength portion (CCD1 with blue collimator) of HIRES 
spectrum of PG1225-079. The number density of absorption lines increases dramatically toward UV wavelengths. 
The absolute flux level and the overall continuum slope are uncalibrated, which is not a problem for the line 
measurements which are always made relative to their local continuum. 
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Fig. 2. — Combined, echelle-order-merged HIRES spectrum of HS 2253+8023 (black), highly smoothed for clarity in 
displaying the broad helium features over the large wavelength range plotted. A T e jy=14,400 K, log g = 8.4 model, 
with element abundances from Tabled! is plotted in red. 
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Fig. 3. — Portion of the HIRES data (black) of HS2253+8023, smoothed by a 5-point boxcar average and displaying 
Ha. Wavelengths are in air and a heliocentric rest frame. A T e jj=14,400 K, log g = 8.4 model, convolved with 
a Gaussian to approximate the instrumental resolution and with the hydrogen abundance from Table 3J is plotted 
in red. The echelle order end occurs in the data near 6553 A, thus the continuum blue- ward of the Ha feature is 
minimally sampled, but is adequate for an EW measurement, which is matched with the model. 
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Fig. 4. — Portion of the HIRES data (black) of HS2253+8023, smoothed with a 5 point boxcar average and displaying 
lines of He I, Ca II, Cr II, and Fe II. Wavelengths are in air and a heliocentric rest frame. A T e j-j=14,400 K, log g 
= 8.4 model, convolved with a Gaussian to approximate the instrumental resolution and with element abundances 
from Table HI is plotted in red. 
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Fig. 7, 
Ni II. 



HS2253+8023. Similar to Figure [U but displaying lines of Ti II, Cr II, Mn II, and a marginal detection of 
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Fig. 8. — Portions of the HIRES spectrum of PG1225-079 are plotted in black, displaying a helium line at He I A5876 
(top panel, data unsmoothed) and the region near He I A4472 including the detection of Mg II A4481 (lower panel, 
data smoothed by a 3 point boxcar average). Wavelengths are in air and a heliocentric rest frame. The upper panel 
spans a region across 2 echelle orders causing the gap near A5896, which is simply filled in with a horizontal line. 
As described in Section 13.21 it is challenging to simultaneously fit the He I lines and achieve an ionization balance 
for iron; the red curve is our best compromise. The downward spike at the He line center is not noise, rather our 
current models do not predict the deep narrow core observed in this star's He I A5876 line, see also Figure [16] and 
the last paragraph of Section [2j The sodium upper limit reported in Table [5] is 0.15 dex higher than that of the 
synthetic Na I D line features near A5890A shown here. Lower Panel: Instead of a broad shallow helium line, a weak 
narrow feature (EW ~10 mA) appears, and is centered at He I A4472 where no high-Z lines are expected. Perhaps 
this is analogous to the narrow core of the He I 5876 line. A weak narrow line also appears in the data at Hel A3889. 
The 10500/8.0 model (green curve) shown for comparison was not used for abundance determinations so the Mg 
abundance of this model is not matched to the data. 
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Fig. 9. — Unsmoothed portion of the combined HIRES spectrum of PG1225-079 is plotted in black, displaying a Mg 
I triplet and many Fe I lines. Wavelengths are in air and a heliocentric rest frame. A T e j^=10,800 K log g = 8.0 
model, convolved with a Gaussian to approximate the instrumental resolution and with element abundances from 
Table is plotted in red. 
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Fig. 10.— PG1225-079. Similar to Figured but for Ti and Cr. 
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Fig. 11. — PG1225-079. Similar to Figure but for V and Ni. The data in the lower panel are smoothed by a 3 
point boxcar average. 
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Fig. 12. — PG1225-079. Similar to Figure EJ but smoothed by a 3-point (upper panel) and 5-point (lower panel) 
boxcar average, and displaying Sc II and Fe I lines; a probable blend of a Ni I feature is also labeled. The apparent 
"mis-fit" of some Fe I lines is due to a combination of the ~ 0.1 dex discrepancy of Fe I and Fe II abundances, 
uncertainties in the modeled line data, and an un-modeled shift of Fe I features (see also Figure [HI) . 



-43- 



0) 



0) 

> 



GO 
0) 



0) 

> 

J— i 




3208 



3213 



3218 



3223 



3228 




3805 3815 3825 3835 3845 3855 3865 3875 
Wavelength (Angstrom) 

Fig. 13. — Unsmoothed portions of the combined HIRES spectrum of G241-6 is plotted in black, displaying Fe II, 
Mg I, Si II and He I lines. Wavelengths are in air and a heliocentric rest frame. Plotted in red is a T e yj=15,300 
K log g = 8.0 model, convolved with a Gaussian to approximate the instrumental resolution, and with element 
abundances from IZuckerman et al.1 (201CJ); the model is blue shifted by -28.0 k m s _1 , the mean velocity from the full 
set of measured absorption lines. See also Figure 4 of IZuckerman et al.1 (|2010l ) for a display of O I lines in G241-6. 
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Fig. 14. — Histograms of (heliocentric) velocity measurements from HIRES data. The velocity bin size is 0.5 km 
s _1 for PG1225 and 1 km s _1 for HS2253; the larger velocity dispersion for the latter star is likely due to the lower 
SNR. Fe I lines are a distinguished set in PG1225, whereas the eight Fe I lines in HS2253+8023 span between -3 km 
s _1 and +12 km s . 
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Fig. 15. — Upper Panel: HIRES observations of Ca II H Sz K emission cores in the white dwarf PG1225-079; zoomed- 
in display of Figure [TJ The local level of the stellar continuum is near 0.55 on this ordinate axis of arbitrary scale. 
The equivalent widths of the core emission features, relative to the lowest depth of the associated calcium lines, are 
~50 mA centered on the peak at A3934.32 A and ~40 mA centered on the peak at A3969.12 A. Lower Panel: Core 
portion of the Ha profile of PG1225-079. All wavelengths are in air and a heliocentric rest frame. 
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Fig. 16. — Narrow cores on the broad He I A5876 line are detected in HIRES spectra of three helium-dominated, 
hydrogen-lined, externally polluted WDs in the ^eff range of 10,000 K to 12,000 K. PG1225-079's observations are 
from Table [TJ GD16 was observed in both February and Nov ember of 2008 with no noted variations of the He line, 
and GD362's observations are from IZuckerman et al.l (|2007l ). The equivalent widths of the core features are: 50 
mA, 45 mA, and 25 mA in PG1225-079, GD16, and GD362, respectively. For this display each WD spectrum was 
Doppler shifted to the laboratory rest frame, by its measured redshift, to show the alignment of the core feature 
with the central wavelength in air A5875.70 A of the helium line. 
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Fig. 17.— Portions of the spectral energy distributions (SED) of PG1225-079 (top) and HS2253+8023 (bottom) 
displaying archival IUE spectra (black, solid line) and GALEX photometry (yellow boxes). The models are convolved 
with 6A Gaussians and scaled to match visual photometry. For PG1225-079 the 10,500 K model is a red, dashed 
line, and the 11,100 K model is a green, dotted line. For HS2253+8023 the 14,000 K model is a red, dashed line, 
and the 14,800 K model is a green, dotted line. An extensive list of UV absorption lines is included in the models, 
with ele ment abundances from the current HIR ES measurements (Tables 0] and [5]) and carbon from previous UV 



analyses (jFriedrich et al 



1999; 



Wolff et al.l 120020 . In both cases the lower temperature model fits the SED well, but 



if there is some unaccounted-for extinction, then the higher temperature model may be appropriate. 
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Fig. 18. — Parent body abundance ratios from Tabled with values f or Earth from lAllegre et al. ( 2001 ) and settling 
times from Tables [4] and 02 50% condensation temperatures are from lLodders ( 2003 ). The fiducial element, Mg, has 
a 50% T c = 1336 K, nearly the same as that of Fe. Oxygen is off the plot at 50% T c = 180 K. The dotted line at 
ordinate value 1 is to guide the eye to Earth-like values. 
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Fig. 19. — Accreted parent body mass compositions of major elements not including hydrogen, for an early phase 
interpretation. For bulk Earth, HS 2253, GD 40, G 241-6 and SDSS 0738, the omission of hydrogen has a negligible 
effect on this display, since it is less than 1% of the bulk mass in each of these systems. Conversely, GD 362 and 
GD 61 have large amounts of hydrogen in their atmospheres, which may or may not be associated with currently 
obser ved high-Z el e ments from polluting parent bodies; s e e Sec tion 16.31 Abu ndance data for th e bulk Earth are 
fr om lAllegre et al.1 (boOlh ; GD362 from IZuckerman et alJ (j2007h : GD4 from klein et al.1 (120101 ): SDSS0738 from 



Dufour et al 



(|2010h : G241-6 from IZuckerman et al.1 (|2010h : GD61 from lFarihi et al.l (j201l|). 
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Fig. 20. — Similar to Figure fl~9l but here in the steady state interpretation of parent body composition. Abundances 
for white dwarf systems are from the same sources as in Figur e 1191 with the element abundance ratios scaled by 
ratios of settling times. Settling time data for GD362 are from iKoesterl (boOflh : HS2253 from Table SJ GD40 from 



Klein et al.l (120101'); SDS S 0738 from our calculations analogous to 
GD61 from lFarihi et all (|201ll ). 



Koesterl feood ): G241-6 are the same as GD40; and 



